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This short article contains the derivation of a general formula for the magnetostatic field in the 
neighbourhood of an axis of symmetry. 



I. INTRODUCTION 

There is a remarkable result in magnetostatic when the field posses cylindrical symmetry, namely, for points 
closed to the axis of symmetry the magnetic field is determined by the distance from the axis and the values it 
takes on this axis. Demonstration of this result appears in many textbooks |lj and is often left as an exercise; 
it consists in giving for each component, the first terms of an ascending series in powers of the distance from 
the axis. 

Whereas the first terms are easily found by successive approximations method, the difficulty rapidly increases 
when we try to derive higher order terms by the same method. The purpose of the following lines is to overcome 
this obstacle by a direct derivation of the general terms of each ascending series. 

II. DERIVATION 

We shall use the cylindrical coordinates (r, 8, z) with the z-axis for the symmetry axis. The magnetic field, 
B, has two components B r and Bg which depend on r and z. It is searched in regions containing the axis of 
symmetry and where current are absent. With this former condition, Maxwell's equations for the magnetic field 
reduce to 

VAB = 0, (1) 
V • B = 0. (2) 

The first equation allows the magnetic field to be derived from a scalar potential 4> defined by 

B = V<j>. (3) 

The substitution of Eq. [3] in Eq. [21 leads to Laplace's equation for the potential 

A<i> = 0. (4) 

Note that the electrostatic field obeys the same equations and the final results for the magnetic field will thus 
be valid for the electrostatic field. 

In cylindrical coordinates the Laplace's equation becomes 

r dr ( dr / dz 2 ^ ^ 

The separation of variables by means of the product functions, 

<Kr,z)=f(r)g(z), (6) 
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leads to the following system of differential equations where k is the separation parameter 

g" + k 2 g = 0, (7) 
r 2f + rf _ r 2 k 2 f = Q (g) 

The integration of the first equation gives e ±lkz . Solutions of the second equation are Bessel's functions of the 
first and second kind with zero order [2j, and whose argument, ikr, is purely imaginary 1 . The Bessel's function 
of the second kind must be rejected because of its singularity at r — 0. We thus retain only Jo, the Bessel's 
function of the first kind 

f(r) = J Q (ikr). (9) 
Now we have the following set of elementary solutions of Eq. (|5j) 



J (ikr)e lkz , (10) 



on which we shall expand the potential 0(r, z). 

This potential is assumed to be a definite function on the z-axis, we write therefore 



(f>(p,z) = <p(z). (11) 

Let's call (f(k), the Fourier transform of <p(z): 

/+oo 
<p(k)e tkz dk. (12) 
-oo 

Since we have Jo(0) = 1, the potential 

/+oo 
<p{k)J {ikr)e ikz dk, (13) 
-oo 

obviously coincides with ip on the z-axis and satisfies Laplace's equation. 

This solution is with no doubt unique. Indeed, if we consider the exterior Dirichlet's problem in which the 
potential is known and equals to ip on the surface of an infinite cylinder surrounding the z axis, the unique 
solution of this problem would tend to our solution when we let the radius of the cylinder tend to zero. 
Along the axis of symmetry, the the magnetic field has the only component B z which is supposed to be known. 
We thus write 

B z (0,z) = b(z). (14) 
The Fourier transform b(k) of b(z) is related to (p(k) by 

b(k)=ik0(k). (15) 
Making use of the expansion of Eq. (|13(l , we derive the component B z 

Bd> f + °° f +OQ 
B z = -^-= I J (ikr)ikif(k)e lkz dk = \ J (ikr)b(k)e lkz dk. (16) 



z 
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Now, we must transform this formula into an expansion containing terms with b(z) and its successive derivatives. 
We have, for this purpose, the following ascending series of Jo 

M*) - E (17) 

n=0 

which we report in Eq. (|16|l after substituting ikr for x; this yields to 

°° C_I r 2Vi p + 00 

B * = Y, 4 |2 / (ik) 2n b(k)e ikz dk. (18) 
The integral is nothing else than the Fourier transform of b^ 2n \z), we thus obtain the desired expansion for B z : 

00 / 1 2\n 

B z = Y,- ZJ ^ Lb(2n) ( z )- (19) 



n! 

n=0 



The -B r component can be derived quite similarly, 



Or 

Since we have J' = —J\ and the following ascending series for J\ 



36 f +oc 

B r = ^f = / ikJ'Skr)v{ky kz dk. (20) 



^W^EVtW' ( 21 ) 

2 n n + 1 

n=o y ' 



we obtain 



00 / 1 2\n 

B z = --Y \~ ir \, fr( 2 " +1 >(z). (22) 
2 ^ n (n + 1) K J y ' 

In particular, the first terms of the expansions of B r and Bg components are easily found 0|: 

B.M-B.M-^)*!^ (23) 

(24) 



BrM) .. (0 ^M + (r!) 
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